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An essential part of experimental program at the future Electron Ion Collider is the study of the 
nuclear structure and dynamics at low x. DIS at low x is characterized by large longitudinal 
coherence length that by far exceeds radii of heavy nuclei. The coherent behavior is essential 
feature of the nuclear matter at low x. This pertains not only to the strong interactions, but also to 
electromagnetic ones. Coherent interactions of a projectile with nucleons in a heavy nucleus are 
characterized by parameters a^A 1 / 3 ~ 1 and aZ ~ 1 in strong and electromagnetic interactions 
respectively. Contributions exhibiting non-trivial dependence on aZ are called the Coulomb cor¬ 
rections. We compute the Coulomb corrections to the cross sections of the semi-inclusive and 
diffractive DIS. We show that they violate the geometric scaling in a wide range of photon vir¬ 
tualities and is weakly ^-independent. In heavy nuclei at low Q 2 the Coulomb correction to the 
total and diffractive cross sections is about 20% and 40% correspondingly. 
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1. Why electromagnetic contribution is important 

DIS at low x is characterized by large longitudinal coherence length that by far exceeds radii of 
heavy nuclei. In QCD, color fields of nucleons in a heavy nucleus fuse to create an intense coherent 
color field which has fundamental theoretical and phenomenological importance. Since nuclear 
force is short-range, only nucleons along the same impact parameter add up to form the coherent 
field. Because the QCD contribution to the scattering amplitude at high energy is approximately 
imaginary, it is proportional to a 2 . Thus, the parameter that characterizes the color-coherent field 
is a 2 A 1 / 3 ~ 1, where A is atomic weight. 

Along with strong color field, heavy-ions also posses strong electromagnetic Coulomb field. 
The electromagnetic force is long range, therefore all Z protons of an ion contribute to the field. 
The QED contribution to the scattering amplitude is approximately real. As a result, the parameter 
that characterizes the coherent electromagnetic field is aZ ~ 1. Since both parameters a 2 A 1 3 and 
aZ are of the same order of magnitude in heavy ions, electromagnetic force must be taken into 
account along with the color one. This observation is a direct consequence of coherence which 
enhances the electromagnetic contribution by a large factor Z. Such electromagnetic contributions 
are called the Coulomb corrections. 

In this article I focus on the Coulomb corrections to the total semi-inclusive and diffractive 
cross sections. At low x the total y*A cross section can be expressed in terms of the total dipole- 
nucleus cross section d as follows (see e.g. [1]) 

Ot/l(x,Q 2 ) = ^j 0 dz j d 2 r<p T/L{r,z)o(x,r), (1.1) 

where Q 2 is the photon virtuality. The light-cone wave functions for transverse and longitudinal 
polarizations of photon are given by 

<$>T = {[r+ (i -z) 2 )a 2 K\(ar)+m 2 f Kl(ar)} , 

f K 

= £ ^4<2V(1 -z) 2 K 2 (ar), 
f K 

where nif is quark mass, z is the fraction of the photon’s light-cone momentum carried by the 
quark, r is the size of the qq dipole in the transverse plane and a 2 = z(l — z)Q 2 + m 2 - 1 The dipole 
cross section d depends non-trivially on A and Z. In the following sections I explain how it can be 
calculated. More details can be found in [2, 3]. 

2. Nuclear dependence of the dipole cross section d. 

At high energies, interaction of the projectile proton with different nucleons is independent 
inasmuch as the nucleons do not overlap in the longitudinal direction. This assumption is tanta¬ 
mount to taking into account only two-body interactions, while neglecting the many-body ones [5], 

1 The relationship between the cross section a = Cj + or and F\ , F 2 structure functions is non-trivial due to large 

Coulomb corrections to the leptonic tensor [4], 


(1.2) 

(1-3) 
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In this approximation we can write the dipole-nucleus elastic scattering amplitude as 

r(b) = i-exp{-Y j {r^(b))}. ( 2 . 1 ) 

where F 11 ( b) is the dipole-nucleon elastic scattering amplitude and (...} denotes average over the 
nucleon position in the nucleus. For the sake of brevity, dependence of the scattering amplitudes on 
x and r is not explicitly indicated. The total dipole cross section can be computed using the optical 
theorem as follows 

<7 = 2 J d 2 blm[ir(b)]. ( 2 . 2 ) 

Strong and electromagnetic contributions decouple in the elastic scattering amplitude at the 
leading order in respective couplings: 

rW=ri 1) +r^ ) . (2.3) 

This is because iT^m is real, while /h!,' is imaginary, as discuss below. Owing to (2.3) we can cast 
(2.1) in the form 

r(fc) = 1 -exp{ -A(r?\b)} -z(r ( ^(b)) } • (2-4) 

Averaging over the nucleus wave function is done as follows 

(ri 1} (*)) =jJ~ dZaj d 2 b a p(b a ,Za)r { P{b-b a ), (2.5) 

where p is the nuclear density. Neglecting the diffusion region, nuclear density is approximately 
constant p = A/(^kR\) for points inside the nucleus and zero otherwise. The range of the nuclear 
force is about a fm, which is much smaller than the radius R,\ of a heavy nucleus. Therefore, h ~ b a 
and 

(ri 1} (b))=j 2y/R* - b 2 xR 2 plf) (0) = H -pT(b ) l - nr 2 a* In ±, (2.6) 

where T{b) = 2 \JR\ — b 2 is the thickness function and p is an infra-red scale. It follows from 

(2.6) thatA(ri 1} ) ~ a 2 A 1 / 3 , which implies that (2.4) sums up terms of order a 2 A 1 / 3 ~ I at a s <A 1. 
Indeed, the leading strong-interaction contribution to the y*N elastic scattering amplitude corre¬ 
sponds to double-gluon exchange. Note also, that the corresponding (;r[* ) is purely imaginary. 

To calculate the electromagnetic contribution to the y'N scattering, note that the proton density 
in the nucleus is Zp/A. Hence 

(r'!,] (*)) = I f—oo dZa f d 2 b a f P (b a ,Za) rl 1 ,) (b - b a ) (2.7) 

= y 2a Jd 2 b a T(b a ) In■ (2-8) 

The leading electromagnetic contribution to elastic y'N scattering amplitude arises from one pho¬ 
ton exchange; the corresponding (/r^j) is purely real. We note, that (2.5) sums up terms of order 
aZ ~ 1 at a <S 1. Had we been interested in purely electromagnetic scattering (e.g. of e e + instead 
of qq) we could have approximated b A> b a ~ R,\ owing to the long-range nature of the Coulomb 
potential. That would have yielded the well-known Bethe-Heitler-Maximon result [6, 7] as shown 
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in [8]. However, in DIS b ^ R A and no such approximation is possible. It follows from (2.2) and 
(2.4) that 

c = lj d 2 b 11 — exp [—A ^r[ 1} (&)^] cos[Z (iT^2 W)]} ■ (2-9) 

Integrals in (2.8) and (2.9) can be analytically calculated in a simple but quite accurate “cylindrical 
nucleus" model (see e.g. [9, 10]), which approximates the nuclear thickness function by the step 
function, viz. T {b) = 2R A if b <R A and zero otherwise. The result is [2] 

a(x,r) = a s (x, r) + 6 em (x,r), (2.10) 

& s (x,r) = 2nR 2 A {1 — exp [~\Q 2 s {x)r 2 ] } , (2.11) 

a em (x, r) = 4^r 2 («Z) 2 ln , (2.12) 

where niM is nucleon mass, W is the y*A center-of-mass energy given by W 2 = Q 2 /x + m 2 N and Q 2 
is the quark saturation momentum (see e.g. [1]). 

Logarithm that appears in (2.12) is the result of integration over the impact parameter from 
R a up to a cutoff b max , which delimits the region of validity of the Weizsacker-Williams approxi¬ 
mation. It is given by b max = max { W 2 z( I — z)/ (»tv( /?1 / + A: 2 )) }, where k is the quark’s transverse 
momentum [8]. The largest size of the qq dipole, corresponding to the smallest k, is ~ I /mj due 
to the confinement. For that reason b max , and hence (2.12), depends on the constituent quark mass 
nif rather than on the much smaller current quark mass m q . 

3. Evolution effects 

Eqs. (2.10)-(2.12) are derived in the quasi-classical approximation where the quark saturation 
momentum Q 2 , and hence the QCD contribution to the total cross section, is x-independent. At 
lower x, such that a s ln( 1 /x) ~ 1, the QCD quantum evolution effects become important and are 
described by the BK equation [11, 12]. According to the BK equation the saturation momentum 
acquires x-dependence in the form Q 2 ~ 4 l/: V U where A is a certain positive number [13]. The 
functional form of the dipole cross section is also evolving with x; (2.11) in that case is the initial 
condition. There are several phenomenological models that describe the evolved dipole cross sec¬ 
tion. In this article I use the Golec-Biernat-Wusthoff model [14] which retains the functional form 
of (2.11) while models the saturation momentum as follows 

Q 2 s=Qo(j)\ (3-D 

where Q a = 1 GeV,x 0 = 3.04• 10~ 4 , A = 0.288, R A = R p A l/3 with R p = 3.1 GeV^ 1 , m f = 140 MeV, 
and Nf = 3. If we neglect the electromagnetic term (2.12) and use (2.11) in (1.1), then we imme¬ 
diately observe that the total y*A cross section exhibits the geometric scaling at Q 2 A> mj. This 
is because x-dependence arises only through the combination r 2 Q 2 (x), and the dipole size r is 
determined by I /Q. 

That the Coulomb correction violates the geometric scaling is evident from (2.12) which, being 
an electromagnetic contribution, does not depend on the strength of the color field determined by 
Q 2 . Unlike the QCD term (2.11), the QED one (2.12) does not evolve much with x. Indeed, at the 
leading-log oder T s / em ~ (1/x) 1+a A™, where A v = 4 In 2{a s N c / k) [15, 16] and A,, m = (I I /?>2)Tia 2 
[17, 18]. Because A em <C A s we can neglect the effect of the QED evolution. 
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4. Total cross section 


Substituting (2.12) into (1.1) and integrating over r we obtain the following analytic expression 
for the Coulomb correction to the total y'A cross section [3] 


n,T /L = (aZ)Hn^^ 8T/L (Q/m f )., 


8t(ti) = 


4 (tj 4 + 7tj 2 + 8) tanh 1 


V \J b 2 +4 

ff +2 


-2t7 V /t 7 2 + 4(?7 2 + 8) r] 3 (t] 2 + 4) 3/2 


gdri) = 4 


t? yf T] 2 +4 (t 7 2 + 6) +4(r7 2 + 3) In 


\/t? 2 +4 


! +4 


T] 3 (r] 2 +4) 


3/2' 


(4.1) 

(4-2) 

(4.3) 


Coulomb correction to the total semi-inclusive cross section (4.1) does not scale with Q 2 JQ 2 and 
therefore explicitly violate the geometric scaling. 

The numerical results are shown in Figs. 1-2. We can see in Fig. 1 that at low Q 2 the QED 
correction for uranium nucleus at x = 10~ 4 can be as large as 20% in the total cross section. It is 
remarkable that the Coulomb correction is non-negligible even at high Q 2 . One should, however, 
interpret cautiously the results of our calculation at high Q 2 since the model that we are using 
does not properly account for the DGLAP evolution. As expected, the relative size of Coulomb 
corrections increases with the nuclear weight and weakly depends on x. 



Figure 1: Ratio of QED and QCD contributions to the total y*A cross section at x = 10 4 as a function of 
Q 2 for silver (solid line), gold (dashed line) and uranium (dotted line) nuclei. 
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Figure 2: Ratio of QED and QCD contributions to the total y*A cross section as a function of Q 1 for gold 
nucleus at x = 10~ 2 (solid line), x = 10~ 3 (dashed line), x = 10~ 4 (dotted line). 
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5. Diffractive cross section 


Total diffractive cross section corresponds to elastic scattering of color dipole on the nucleus. 
It can be written as 


T diff 
’T/L 


yX 'Q 2 ) = In h dl I 6 el (x,r), 


where the total elastic dipole-nucleus cross section reads 

d e, (x , r) = Jd 2 b 11 - exp [—A (iF s ) - Z (iT em )]\ 2 
A simple calculation yields [2] 

6 eI (x, r) = of (x, r) + 6 em (x, r ), 

where o em is the QED contribution given by (2.12), while the QCD contribution is 

1 


of (x. r) = nR\ l 1 — exp 




(5.1) 


(5.2) 


(5.3) 


(5.4) 


Numerical results are shown in Figs. 3-4. We can see in Fig. 4 that at low Q~ the QED 
correction for uranium nucleus at x = 10 4 can be as large as 40% in the diffractive cross section. 
The Coulomb corrections to the diffractive cross section are much larger than in semi-inclusive 
one, see Fig. 3. The reason is that at the leading order in coupling, electromagnetic interaction is 
elastic. 
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Figure 3: Ratio of QED and QCD contributions to the total (solid line) and diffractive (dashed) y' ,4 cross 
section as a function of x for gold nucleus at Q 2 = 1 GeV 2 . 


6. Conclusion: beware of Coulomb corrections 

An essential part of experimental program at the future Electron Ion Collider (EIC) is the study 
of the nuclear structure and dynamics at low x. Nuclear matter at low x exhibits highly coherent 
behavior, the most striking manifestation of which is the geometric scaling [19, 20, 21, 22, 23]. 
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Figure 4: Ratio of QED and QCD contributions to the diffractive y'A cross section at x — 10 4 as a function 
of Q 2 for silver (solid line), gold (dashed line) and uranium (dotted line). 


In semi-inclusive DIS the geometric scaling means that the total y *P and y *A cross sections scale 
with a dimensionless ratio Q 2 /Q 2 (x), where Q 2 is photon virtuality and Q s (x) is the quark satu¬ 
ration momentum. While the geometric scaling is a signature of coherence in strong interactions. 
Coulomb corrections is a manifestation of coherence in electromagnetic interactions. Moreover, 
in DIS off heavy nuclei the two coherence effects are delicately entangled and separate only in a 
crude cylindrical nucleus model. 

Results presented in this work indicate that Coulomb corrections play an important role in 
the low x DIS off heavy nuclei in a very wide range of Q 2 and x. In particular, they violate the 
geometric scaling. Since the focus of the EIC program is on the strong interactions, one would 
like to learn how to subtract the Coulomb corrections from the measured quantities. We showed 
that Coulomb corrections can in principle be reliably computed. Since we used an oversimplified 
model for nuclear wave function that allowed us to derive analytical formulas, it is important to 
refine our estimates by using more realistic nuclear profiles and more sophisticated low x evolution 
models. 
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